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Absorption spectra of closed-shell Na2 , Na J , Na4 , Na^" , Na6 , Na^ , and Nas clusters are calculated 
using a recently implemented conserving linear response method. In the framework of a quasiparticle 
approach, we determine electron-hole correlations in the presence of an external field. The calculated 
results are in excellent agreement with experimental spectra, and some possible cluster geometries 
that occur in experiments are analyzed. The position and the broadening of the resonances in 
the spectra arise from a consistent treatment of the scattering and dephasing contributions in the 
linear response calculation. Comparison between the experimental and the theoretical results yields 
information about the cluster geometry, which is not accessible experimentally. 

PACS numbers: 78.67.-n Optical properties of low-dimensional, mesoscopic, and nanoscale materials and 
structures; 36.40.Vz Optical properties of clusters; 73.22.-f Electronic structure of nanoscale materials and 
related systems; 36.40. Cg Electronic and magnetic properties of clusters 



I. INTRODUCTION 

The electronic configuration of small metal clusters is 
closely connected with their geometrical, chemical and 
optical properties^. Different aspects of this interrelation 
for small metal clusters have been extensively studied ex- 
perimentally 2 — and theoretically-^- —. For instance, the 
measured photoabsorption cross section of mass-selected 
clusters yields important information about the inter- 
play between the underlying nuclear configuration and 
the electronic structure. However, the geometry of the 
clusters in the gas phase is not directly accessible in ex- 
periments. Ab initio theory can help to determine the 
structures by choosing, among many possibilities, the 
most stable geometrical configuration, for which the cal- 
culated absorption spectrum is in good agreement with 
the experimental results. 

Early theoretical attempts to understand the nature 
of electronic correlation in small alkali-metal clusters 
are based on the shelU 3 - and density functional theory 
(DFT)-based jelliumi^"— models. Quantum molecular 
approaches based on the Configuration Interaction (CI) 
procedures are powerful theoretical tools capable of yield- 
ing accurate electronic excited states and important in- 
formation about the nature of the optically allowed tran- 
sitions between them. The multi-reference single and 
double excitation CI has been successfully used to ac- 
count for the dominant features of the absorption spec- 
tra of small sodium clusters^"—, calculated from the os- 
cillator strength of the transitions between the ground 
and the excited states. Also, time-dependent (TD) DFT 
methods are widely employed to compute the photoab- 
sorption cross sections of small metal clusters, with var- 
ious levels of sophistication for the exchange-correlation 
functional 2 ^ - —. A different way to compute optical spec- 
tra is to make use of the quasiparticle picture, where 



a typical photoabsorption process involves an incoming 
photon and the transition of an electron from an occu- 
pied to an empty quasiparticle state. The determination 
of the cross section can be achieved by a calculation of a 
generalized "four-point" (i.e., two-particle) electron- hole 
correlation function, which obeys a Bethe-Salpeter equa- 
tion (BSE)22r— . In practice, this is usually done in a 
two-step calculation: first the DFT or Hartree-Fock (HF) 
single-particle states are used as input into a GW-like 
procedure to obtain the quasiparticle corrections to the 
Kohn-Sham or HF ground- state eigenvalue s 35 ^ 6 . Then a 
two-particle BS effective equation with a screened direct 
and an unscreened exchange electron-hole interaction is 
solved. While the first step properly describes charged 
excitations (i.e., electron removal or addition, specific to 
direct or inverse photoelectron experiments), the second 
step accounts for neutral excitations (i.e., electron- hole 
pair creation, specific to photoabsorption experiments). 
A comprehensive comparison between the TDDFT and 
the GW+BSE approaches can be found in Ref. [32|. To 
facilitate the numerical procedure, the screened Coulomb 
interaction, which is obtained in the GW step and which 
enters the BSE, is most often taken to be the statically 
screened Coulomb interaction. This simplification allows 
one to cast the effective two-particle BSE in the form of 
an energy-independent eigenvalue problem, which yields 
the excited states corresponding to neutral excitations of 
the system. As in the case of CI, the optical absorp- 
tion cross section is then determined by the oscillator 
strengths of the transitions to the excited states, and 
a phenomenological parameter for the resonance peak 
broadening is added in order to compare the calculated 
spectra with the shapes of the experimental peaks. This 
phenomenological time for the individual transition mim- 
ics the combination of thermal broadening and electron 
lifetime. 
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We have developed, in a quasiparticle picture, an 
alternative method to determine the photoabsorption 
spectra of small sodium clusters. We employ a first- 
principles based linear response calculation for the dy- 
namical electron-hole correlation function in the pres- 
ence of an external potential. Although we do assume 
a "background" quasiparticle lifetime, the finite width 
of the absorption peaks contains electron-hole correla- 
tion contributions beyond a mere convolution of single- 
particle lifetimes. A novel aspect of our approach is that 
we use a quantum kinetic equation for the electron-hole 
correlations driven by an external coherent optical field 
to obtain a closed equation for the two-particle BSE with- 
out the need of a four-point correlation function. Techni- 
cally, we include mean-field HF contributions (accounting 
for direct and exchange bare Coulomb interaction) to- 
gether with correlation contributions beyond mean-field, 
described by complicated integral kernels. An impor- 
tant difference from the two-step GW+BSE approaches 
mentioned above is that we do not include correlations 
on the single-particle level and then solve the BSE with 
an effective two-particle interaction described in terms 
of a simple integral kernel. Instead, we follow the orig- 
inal prescription of deriving conserving approximations 
for transport equations by Baym and Kadanoff. We thus 
obtain a BSE for the electron-hole correlation function 
that is consistent with the electronic single-particle prop- 
erties^. This type of conserving approach (in the sense of 
Baym and Kadanoff) avoids possible problems of double- 
counting correlation contributions which may arise in a 
two-step GW+BSE-like calculation due to the separate 
treatment of single- and two-particle correlations. As 
already shown by Baym and Kadanoff, such a conserv- 
ing approach fulfills important sum-rules by construction. 
For instance, the particle- number conservation law for 
the electron-hole correlation function, which can be read 
as a sum rule for the absorption cross section in finite 
systems, serves as a check of our numerical calculations. 



II. THEORY 

The BS approach employed to calculate the absorption 
spectra is described in detail in Refs. |38| and [39|, and is 
briefly reviewed in the following. The central quantity of 
this approach is the retarded density-density correlation 
function x r ( r ? t] r', £'), which describes the linear response 
of the system to an external potential U(r,t). For finite 
systems, we use the matrix element 



(nin 2 |x r (^^)l n 3^4) 



SU n3n4 (t') 



(1) 



u=o 



in a basis of HF eigenfunctions {ip n (r)}. Here, n labels 
the HF spin orbital, and p n2ni (t) = (c^ 2 (i)c ni (t)) with 
the annihilation (creation) operators c n (cjj for an elec- 
tron in HF orbital <p n . In linear response, the electronic 
distribution functions p nn are unchanged, and the only 



nonzero elements in Eq. ([T]) are those with indices that 
pair occupied and unoccupied states, i.e., Eq. (pQ) is the 
electron-hole correlation function. 

The photoabsorption cross section is calculated from 
the Fourier transform of the imaginary part of x r ac- 
cording to 



e c 



^2 Im(n 2 ni|x r M|n 3 n 4 ) 



(2) 

where d nin2 = J y?* (er) cp n2 d 3 r are the electric dipole 
matrix elements in the molecular orbital representation 
and e is the electron charge. 

The two-particle BS equation for x r is 
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HF 



)(^1^2|X r (^)|^3^4)- 
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(nin2|v|n 5 n 6 ) J (n 5 n 6 |x r '(^)|n 3 n 4 ) 

^ (n 1 n 2 \A(u)\n 5 n 6 )(n 5 n 6 \x r (uj)\nsn4) = (3) 



n 5 n 6 

with the correlation kernel 



(nin 2 |A(cj)|n 3 n 4 ) = 

,/ni(l — fn 5 )fn 6 0- ~ «/ni)«/n 5 (l — fn 6 ) 



x (nin 5 |v|n 3 n 6 ) (n 6 n 2 \v\n 5 ri4} 
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/n 7 (l — fn 6 )fn 5 + (1 — fn 7 )fn 6 {^- ~ fn 5 ) 



x {ri4n 5 \v\n 7 n 6 ) (n 6 n 2 \v\n 5 n 7 ) 
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In Eqs. © and @j) f n = p nn are the level occupation 
numbers assumed to be equilibrium Fermi-Dirac distri- 
butions at T = OK. The bare Coulomb matrix element 
is defined by 

(nin 2 |v|n 3 n4> = / d 3 nd 3 r 2 ^ (ri)y£ 2 (r 2 ) (5) 
x v(|n -r 2 |)^ n3 (ri)^ n4 (r 2 ). 

The denominators in the expression for the correlation 
kernel (|4]) contain the quasiparticle energies including an 
inverse lifetime e n = s n + ij n . In principle, one can 
use for e n the results of a quasiparticle calculation such 
as a GW calculatio n 38 ! 40 ! 41 , but here we employ the HF 
energies and a constant broadening for all states whose 
magnitude is in agreement with GW calculations. We 
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stress that introducing a fixed quasiparticle lifetime for 
the calculation of A (a;) does not mean that the broaden- 
ing of the resonances in x r is given by this value. Rather, 
the resonance broadening, which is responsible for the fi- 
nite width of the peaks in the absorption spectrum, is 
due to the imaginary part of the full correlation contri- 
bution A((j). The real part of A (a;) yields the shift in 
the resonance energies. Both effects, renormalization and 
broadening of the x v ( (j0 ) spectrum, are thus related. 

Equation [3] for the electron-hole correlation can be 
solved at different levels of approximation. Neglecting 
all Coulomb interactions and setting 



(nin 2 |A(cj)|n 3 n 4 ) 



(6) 



one obtains for x r the independent particle (Lindhard) 
result 



(nm 2 \xl(u)\nzn A ) 



h(u + irf) 



„HF 



(?) 



Including direct Coulomb matrix elements, exchange ele- 
ments and correlation contributions in that order, Eq. (|3]) 
resembles a BSE in the ladder approximation (denoted 
by BS-L in the following) plus exchange (BS-LX) and 
correlation terms (BS-LXC), respectively 38 . Recall that 
we use loosely here the term BSE for Eq. ([3]), even though 
the BSE is usually written for the four-point generaliza- 
tion of x- 



III. PHOTOABSORPTION SPECTRA OF 
SMALL NA CLUSTERS 

In this section we discuss photoabsorption spectra of 
closed-shell Na 2 , Na^, Na 4 , NaJ, Na6, Na^ , and Nas 
clusters. The numerical results are calculated accord- 
ing to Eq. ([2]) and compared with experiment. We 
choose clusters with closed-shell configurations with an 
even number of electrons, for which all molecular orbit als 
are doubly occupied or empty. Moreover, Na 2 , Na^, 
and Nag are magic- number clusters, for which the num- 
ber of valence electrons equals the spherical-shell closing- 
numbers 2 and 8, respectively 6 . The electron- hole corre- 
lation function is determined by solving Eq. (|3]) with the 
correlation kernel (|4j). 

To obtain the equilibrium (T = OK) geometrical 
configuration, we perform a structural optimization us- 
ing the HF approximation, which yields the single- 
particle ground state energies and wave functions. For 
the Na atoms we use the lanl2dz basis set (double- 
zeta set with relativistic Los Alamos effective core po- 
tential) 42 , in which the 3s 1 valence electron of each 
Na atom is represented by contracted Gaussian-like 
atomic orbitals (3s3p/2s2p) while the contribution of 
the core electrons is treated using effective core poten- 
tials (ECP). To test the stability of our results with re- 
spect to the basis set used in the calculations, we have 
also performed calculations with six other basis sets: 



SHC (Goddard/Smedley ECP) 43 , CEP-4G (minimal 
set and Stephens/Basch/Krauss ECP), CEP-31G (split 
valence set) and CEP-121G (triple-split basis set) 44 , 
LANL2MB (MBS set and Los- Alamos ECP) 42 , and SD- 
Dall (double-zeta set and Stuttgart/Dresden ECP) 45 . 
For larger basis sets such as lanl2dz, SHC, CEP-121G 
and SDDall the Na-Na distances vary by less than 1% 
at the end of the geometry optimization, and the peak 
positions of the absorption spectra differ by less than 
0.1 eV. For smaller basis sets the Na-Na distances vary 
at most by 0.3 A while the peak positions have devi- 
ations of less than 0.3 eV. Overall, a good comparison 
with experiment is found for lanl2dz, and in the fol- 
lowing we refer only to results obtained using this basis 
set. When the absorption spectra are calculated from 
the density-density correlation function obtained in the 
Lindhard, BS-L and BS-LX approximations, we use a 
constant broadening T = 0.03 eV in Eq. (|6]). When corre- 
lation terms A are included in the calculations we use for 
the quasiparticle inverse lifetime the value 7 = 0.15 eV 
in Eq. gj). 
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FIG. 1. Experimental and calculated absorption spectra of 
Na2. (a) Experimental spectrum adapted from Ref. Kg. (b) 
Cross section calculated from % r m the BS-LXC (solid line) 
and BS-LX (dot-dashed) approximation. The inset shows the 
cluster structure with the distance in A. 

Within the present approach, the computed spectra 
should fulfil the /-sum rule for finite systems 

x Im(n 2 ni|x r M|n 3 n4) ~^N e 

for the conserving BS-LX and BS-LXC approximations, 
Here, m is the electron mass and N e is the total number 
of electrons in the system. Equation JH} is equivalent 
to the Thomas- Reiche-Kuhn sum rule for the absorption 
cross section and is directly related to the particle number 
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conservation law2£. For our results, Eq. (jSj) is numerically 
checked and the /-sum rule is fulfilled to better than 99%. 

Calculated photoabsorption spectra for Na2 are shown 
in Fig. [T] and compared to experimental data. Only 
when correlations are included in the computation of 
X r , the theoretical spectrum reproduces the experimen- 
tal absorption peaks which occur around 1.97, 2.55, and 
3.75 eV. The pronounced difference between the BS-LX 
and the BS-LXC results (for the same ground-state clus- 
ter geometry which is optimized) shows the importance 
of electronic correlation contributions to the absorption 
spectra. 
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FIG. 2. Experimental and calculated absorption spectra of 
NaJ . (a) Experimental spectrum adapted from Ref. |9|. (b) 
Cross section is calculated in the BS-LXC (solid line) and BS- 
LX (dot-dashed) approximation. The inset shows the cluster 
structure with the distances in A. 



for Na^". Again, the theoretical spectrum resolves both 
the larger experimental peak at 2.62 eV as well as the 
smaller peak at 3.33 eV only when correlations are taken 
into account. Without correlation contributions, the BS- 
LX approximation incorrectly describes the positions of 
the peaks. 

In the case of Na4, the BS-LX and BS-LXC approxi- 
mations yield peaks in the absorption cross section that 
are energetically close, as can be seen from the panel 
(b) of Fig. [3j However, the inclusion of correlation ef- 
fects results in a better agreement with the experimen- 
tal spectrum shown in panel (a) of Fig. [3l for which the 
main peaks are centered around 1.8 and 2.5 eV. Panel (c) 
shows the calculated spectra corresponding to the Lind- 
hard and the BS-L approximation for x r - These approxi- 
mations yield completely different peak positions that are 
red shifted with about 2 eV. Moreover, the /-sum rule is 
not fulfilled for the Lindhard and the BS-L approxima- 
tions, due to the inconsistency between the single- and 
the two-particle quantities used. 
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FIG. 3. Experimental and calculated absorption spectra of 
Na4. (a) Experimental spectrum adapted from Refs.0 and0. 
(b) The cross section obtained from the BS-LXC (solid line) 
and BS-LX (dot-dashed) approximation. The inset shows the 
cluster structure with the distances in A. (c) The spectra ob- 
tained from the Lindhard (dashed) and BS-L (dotted) ap- 
proximations. 
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FIG. 4. Experimental and calculated absorption spectra of 
Na^". (a) Experimental spectrum adapted from Refs. and 
Hq . (b), (c) Cross sections in the BS-LXC approximation 
for clusters with different symmetries: C2v (solid line), C4v 
(dotted), D2d (dot-dashed), and (dashed). The insets 
show the cluster structures with the distances in A. 



To investigate possible geometries of the clusters that 
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TABLE I. Ground-state total energies (in eV) of the Na^, 
Na6, Na^, and Nas clusters for several optimized geometrical 
structures (shown in the panels (b) and (c) of Figs. SI El El 
andQ respectively). The configuration with the lowest energy 
corresponds to the most stable geometry. 

Najj" Na 6 Na^ Na 8 

-21.605 (D 2d ) -29.970 (planar) -31.866 (D 5h ) -40.121 (D 2 c d } ) 
-21.476 (C 2v ) -29.821 (C 5v ) -31.329 (C^ } ) -39.934 (D 4d ) 
-21.168 (D 3h ) -29.526 (D 4h ) -30.984 (C 3v ) -39.661 (D^ } ) 
-20.526 (C 4v ) -28.937 (C 2v ) -30.829 (C 2 C V } ) -39.934 (D 4d ) 



may play a role under typical experimental conditions, 
we consider four different symmetries for Na^~: C2 V7 CW, 
D 2 d, and D 3 h- In Fig. 2] we present spectra calculated in 
the BS-LXC approximation for each symmetry, and com- 
pare them to experiment. As can be seen from the panels 
(b) and (c) of Fig.[4j the cluster structures which yield the 
best agreement to the measured cross section correspond 
to the D2d and C2 V symmetries. They best resolve the 
experimental peaks at 2.2, 2.7, and 3.3 eV. The D2d and 
C2v have similar total ground state energies (0.13 eV dif- 
ference), and are energetically more stable than the 
and C4 V structures. In an experiment, at finite tempera- 
ture, it is possible that both the C2 V and the D2d clusters 
are responsible for the shape of the individual peaks, but 
a clear distinction between the weight of the contribution 
of each cluster is not possible with a T = K theory. 
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FIG. 5. Experimental and calculated absorption spectra of 
Na6- (a) Experimental spectrum adapted from Ref. 47. (b), 
(c) Cross sections in the BS-LXC approximation for clusters 
with different symmetries: planar (solid line), Cs v (dotted), 
C 2v (dot-dashed), and D 4 h (dashed). The insets show the 
cluster structures with the distances in A. 

For Na6 we have also investigated four different sym- 
metries of the cluster structure: planar, Cs v , C2 V , and 
D^. From Fig. [5] one can see that the cross section 



of the planar and the C2 V geometries agree best with 
the experimental spectrum. The energetically most sta- 
ble structure is the planar Na6- The planar structure 
overestimates by 0.1 eV the position of the dominant ex- 
perimental peak at 2.1 eV. It resolves the peak around 
2.7 eV, but predicts another peak at 1.75 eV, which is 
not present in the experimental spectrum. In the case 
of the C2v cluster, the most prominent peak agrees well 
with the experimental peak at 2.1 eV, while the posi- 
tion of the smaller experimental peak around 2.7 eV is 
underestimated by the calculation by 0.2 eV. 
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FIG. 6. Experimental and calculated absorption spectra of 
Na+ . (a) Experimental spectrum adapted from Ref. flOl. (b), 
(c) Cross sections in the BS-LXC approximation for clusters 
with different symmetries: two C 2v (solid and dot-dashed 
lines), Dsh (dotted), and C3v (dashed). The insets show the 
cluster geometries with the distances in A. 

In Fig. [6] we show a comparison between the exper- 
imental and the theoretical absorption spectra of Na^ , 
calculated for two different clusters of C2 V symmetries 
and also for the Dsh and C3 V geometries. The struc- 
ture with the most stable energetic configuration is the 
Dsh, which also reproduces well the 2.4 and 3.3 eV ex- 
perimental peaks, although with a small offset of 0.1 eV. 
The experimental peak at 2.4 eV is well reproduced by 
the structures from panel (c) of Fig. [6l while the C2 V 
structure from panel (b) resolves well the experimental 
peak at 2.65 eV. 

The measured spectrum of the Nag cluster in Fig. [3 ex- 
hibits a single broad peak extending from 2.2 to 2.8 eV. 
The theoretical result is computed for four different struc- 
tural symmetries: D4d, D2h, and two different D2d struc- 
tures, as shown in Fig. Both the D4d and the D2d 
structures from panel (b) give a pronounced peak around 
2.55 eV, which agrees well with the experimental data. 
However, the D2d cluster predicts a smaller peak at 
2.19 eV not present in the experiment. The D2d (which 
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FIG. 7. (color online) Experimental and calculated absorp- 
tion spectra of Nas- (a) Experimental spectrum adapted from 
Refs.0and0. (b), (c) Cross sections in the BS-LXC approxi- 
mation for clusters with different symmetries: D4d (solid line), 
two D 2 d (dotted and dot-dashed) and D 2 h (dashed). The in- 
sets show the cluster structures with the distances in A. 



has the lowest ground state energy) and the D2h struc- 
tures from panel (c) yield several photoabsorption reso- 
nance peaks that lie in the same energy interval as the 
measured spectra. Therefore, different possible configu- 
rations of the cluster geometries and finite temperature 
might be responsible for the strongly broadened experi- 
mental line. However, as mentioned above, it is not pos- 
sible to assign a clear weight for each individual structure 
computed for T = OK. At finite temperatures phonons 
are present and the eigenmodes of these vibronic states 
do not only broaden the peaks, but they may also alter 
the selection rules, leading to new peaks in the absorption 
spectr a 48 ! 49 . 

Summarizing, the calculated BS-LXC spectra pre- 
sented here are in very good agreement with the mea- 
sured photoabsorption cross sections, given the potential 
problems for experiment-theory comparisons which are 



due to the unknown cluster structures, finite tempera- 
ture and inhomogeneous broadening effects that tend to 
wash out the intrinsic spectral features of the resonances. 



IV. CONCLUSIONS 

The absorption spectra of small closed-shell Na clusters 
were calculated by means of a linear response approach 
for the electron-hole correlation function, defined as the 
functional derivative with respect to a weak external 
perturbation of the non-equilibrium single-particle den- 
sity matrix. We numerically solved a Bethe-Salpeter-like 
equation for electron-hole correlation that obeys physical 
conservation laws by construction. In particular, we nu- 
merically checked that the /-sum rule is fulfilled to better 
that 99%. We found that correlations effects beyond the 
mean- field limit need to be included in the calculations in 
order to properly account for the positions and the widths 
of the resonance peaks. In the present approach, the fi- 
nite broadening of the photoabsorption lines is attributed 
to the correlations between the interacting quasielectrons 
and holes with a finite lifetime. 

The computed cross sections of the clusters displayed 
an overall good agreement with the measured spectra 
with respect to the positions and the shapes of the pho- 
toabsorption lines. For the smaller clusters Na2-Na3 we 
found the geometrical configurations that yield theoreti- 
cal spectra in excellent agreement with the experiment. 
For the larger clusters Na^-Nag we have considered dif- 
ferent energetically stable geometries, in order to ac- 
count, within a T = OK theory, for finite-temperature 
effects which are unavoidable in experiments. We were 
able to pinpoint the cluster geometries likely present in 
the experiments, where the structure cannot be directly 
obtained. 
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